Transport through correlated quantum dots: An investigation using the functional 

renormalization group 
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Calculations using the (exact) fermionic functional renormalization group are usually truncated 
at the second order of the corresponding hierarchy of coupled ordinary differential equations. We 
present a method for the systematic determination of higher order vertex functions. This method 
is applied to a study of transport properties of various correlated quantum dot systems. It is shown 
that for large Coulomb correlations higher order vertex functions cannot be neglected, and a static 
approximation is insufficient. 
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I. INTRODUCTION 

Recently, electron transport through ultrasmall 
fermionic systems (quantum dots) has been investigated 
intensivelji^i^i^i^. Motivation for these studies is the po- 
tential application of quantum dots to quantum comput- 
ers and quantum measurement devices^. The transport 
process can be described using the concept of Coulomb 
blockade, when Coulomb correlations are large com- 
pared to other energy scales in the quantum dot sys- 
tem. However, a detailed analysis of the conductance and 
level occupancies of the quantum dot states as a func- 
tion of an applied gate voltage shows a rich structure, 
which is related to Kondo physics. The theoretical de- 
scription of such effects requires non-perturbative meth- 
ods, and many studies employ Wilson's numerical renor- 
malization group (NRG). Recently, functional renormal- 
ization group techniques have emerged as a new non- 
perturbative tool to describe mesoscopic phenomena the- 
oreticallySii. 

While the functional renormalization group equation 
(fRG) is exact in principle, for practical applications it 
must be truncated^. In particular, most fRG investi- 
gations of mesoscopic systems use a static approxima- 
tion (i.e. energy independent self-energies and higher 
vertex functions), and the hierarchy of differential equa- 
tions corresponding to the fRG is truncated at second 
order. Since fRG calculations truncated at second or- 
der are quite 'cheap' numerically in comparison to NRG 
or DMRG (density matrix renormalization group) cal- 
culations, it is possible to map out the potentially huge 
parameter space of quantum dot systems with reasonable 
effort. Functional renormalization group calculations be- 
yond the static approximation are technically difficult 
and require a significant numerical effort as is demon- 
strated in model studies for single quantum particles^. 

The uncertainties incurred by the truncation of the 
fRG are not easily controlled. Therefore, it is worth 
studying the static approximation beyond second order 
truncation. Moreover, methods to go beyond the static 
approximation with reasonable effort should be devel- 
oped, and the limits of the static approximation should 



be assessed. 

It is the purpose of the present paper to investigate 
the influence of vertex functions beyond second order in 
a static approximation. To this end, we develop a method 
which enables to extract the complete set of ordinary dif- 
ferential equations from the underlying functional differ- 
ential equation (section II). Our method uses a different 
regularization scheme than the one usually employed for 
fermionic systems, which reg ularizes the free propaga- 
tor (see e.g. Refs. |^ and Il0[ ). Our cutoff procedure is 
adapted from the scheme used in Ref. 0, which adds a 
regulator to the action. We are able to show that the 
two procedures yield the same set of ordinary coupled 
differential equations for the standard hard cutoff regu- 
lator. Furthermore, it is shown that this set is finite in 
principle for a fermionic system. (This is in contrast to 
Bose systems where this set is infinite.) The number of 
differential equations obtained nevertheless grows expo- 
nentially with the number of degrees of freedom of the 
mesoscopic system under consideration. 

For quantum dot systems described by only a small 
number of electronic states, a study considering all pos- 
sible vertex functions is feasible (within the constraints of 
the static approximation). In this way, the limits of the 
static approximation can be assessed. In order to do that 
we investigate transport properties of various quantum 
dot devices with different couplings to the external leads 
(section III). It is shown that consideration of higher or- 
der vertex functions is important in order to describe 
transport properties of correlated quantum dots quanti- 
tatively, in particular for systems with strong Coulomb 
correlations. However, in many cases we find that the set 
of differential equations to be solved develops singulari- 
ties if we increase Coulomb correlations beyond some crit- 
ical value, and the renormalization flow becomes unphys- 
ical. This is, of course, not unexpected, and just shows 
that for a quantitative analysis of mesoscopic systems 
with strong Coulomb correlations one eventually needs 
to go beyond the static approximation. This will be ad- 
dressed in a forthcoming publication using the methods 
developed here but including a wave function renormal- 
ization in the kinetic energy term of the action. 
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II. FUNCTIONAL RG SCHEME FOR 
INTERACTING FERMIONS 



The effective average action for interacting 

Fermions evolves according to the functional renormal- 
ization group equatio: 



7.11 
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dk 
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dk 



(1) 



Here, we assume that the space of possible states of the 
Fermions is suitably discretized, so that a A'^-component 
vector of Grassmann variables 0(t) = ((/)i(r), . . . , (I)m{t)) 
describes the evolution of the interacting Fermion system 
in imaginary time t. Each index a G {!,..., iV} repre- 
sents all quantum numbers necessary to completely spec- 
ify a Fermionic state, e.g. spin projection a and position 
j in a linear chain of electrons. The functional deriva- 

(2) 

five of Ffc with respect to </>* and </> is denoted by F^ \ 
The regulator Rk is introduced in order to suppress ther- 
mal and quantum fluctuations at energy or momentum 
scales k larger than any physical scale relevant for the 
problem that is beeing investigated. With decreasing fc, 
the regulator gradually "switches on" such fluctuations 
until they are fully included at fc = 0, i.e. at fc = the 
regulator Rk vanishes. A concrete choice for Rk will be 
discussed below. The initial condition for the evolution 
described by Eq. ([T]) is obtained from the Hamiltonian 
defining the problem to be solved^. The trace in Eq. ([T|) 
is to be performed over all relevant quantum numbers. 

In order to solve a functional differential equation like 
Eq. ID) in practice, it must be truncated. This essen- 
tially entails a suitable transformation of the functional 
differential equation into an infinite set of coupled partial 
or ordinary differential equations. This set of differential 
equations is then truncated at a suitable order. The stan- 
dard technique expands both sides of the equation for the 
effective average into a Taylor series about a r indepen- 
dent vector 00 using 0(t) — 0o + '*/'(''")• One then obtains 
an infinite set of ordinary or partial differential equations 
for the Taylor coefficients of the various powers in ij^ir). 

For most practical applications, the truncation proce- 
dure described so far is still too general, and we need 
further simplifications: One standard way to proceed is 
to prescribe the functional form of the effective average 
action more specifically: In this paper we assume that 
the effective average action takes the form 



Tk\ 



At 



N 



d 



a— 1 

(2) 

where the 'effective potential' Uk does not depend on 
derivatives of the Grassmann variables with respect to 
the imaginary time r. Assuming the specific form 
dH) yields energy independent (static) self-energies and 
higher order vertex functions. Of course, more gen- 
eral functional forms for F^ lead to energy dependent 
vertex functions [e.g. a wave function renormalization 



Zk {4>* (t), (f){T)) multiplying the kinetic energy term in 
Eq. but such a truncation will not be investigated 
in this paper. 

Since the effective potential Uk must be a Grassmann 
scalar, its form in terms of </i(t) is fairly well determined: 
It must be a linear combination of all possible products of 
elements of the set of Grassmann variables {0^, 0q} with 
an even number of factors. Due to the nil-potency of the 
Grassmann variables the number of possible products is 
finite, i.e. the number of terms in the linear combinations 
is finite. Since the number of particles is conserved in the 
physical models that we want to investigate, we only need 
to include terms with an equal number of 0* and (pa, 
which correspond to Fermion creation and destruction, 
respectively. 

Here is an example for N — 2: Due to the nil-potency 
of the Grassmann variables, the effective potential has 
the general form 

Uk{<j}*,4>) = flo.fc + '/'I'/'iaii.fe -I- 02</'2a22,fe + (3) 

(/il 020X2, fc + 0201021, fe + 0l'/'2'/'l'/'2ai212,fe- 

The fc-dependent coefficients in this expression will be 
called 'running couplings'. The running couplings have 
direct physical significance: Since the effective average 
action at fc = is the generator of the vertex functions, 
each running coupling corresponds to an n-point vertex 
function, where n corresponds to its number of Grass- 
mann factors. In particular, oq is directly related to the 
ground state energy of the Fermion system, and the aij^k 
correspond to the self-energies. For simplicity, the k- 
dependence of the couplings will not always be explicitly 
indicated. 

Following the general procedure outlined above, we 
will now determine flow equations for the running cou- 
plings. To this end we first expand the effective aver- 
age action ^ about r-independent Grassmann variables 
00 = (001, ■ • ■ , 0ow), i.e. 4>a[T) = 00a + V'a(T) with 
0Oq = 0Q (0) . Up to second order in ipa (t) , the expansion 
is given by 

r,[0*,0] = /3t/f)(0S,0o) 

1 rP rP 



+ ^ dr dr'*t(r)7f (05,0o)*(r') 
+ ... 



(4) 



with*(T) = (^i(r),...,Vjv(T),V^i*(r),...^Ar*(T)). The 
lowest order term is determined by the effective potential 

f/r(0S,0o) = t/fe(0*(r),0(r))|^.^^^o. (5) 

The second order term contains the second derivative of 
the effective potential 



M = (E'9, + C/f' (0S,0o))<5(r-T'), 



d'Uk 



(6) 
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with E' = diag(l, . . . , 1, —1, . . . , —1). The term in first 
order in ^! drops out. For example, for N = 2 assuming 



the symmetry ai2 
the matrix 



a2i, one easily finds from Eq. Q 



J 



(2) 



/ ail - 0(5200201212 Ol2 
ai2 + 00100201212 022 

^ 

V -05x00201212 



'01 



0O20O1O1212 
00100101212 
''O2O12I2 







"0O10O2O1212 
-On + 0^2 00201212 
— O12 — 0o20OlOl212 



-O12 
-O22 



00100201212 



^0100201212 
^^0100101212 



(7) 



r 



In order to do the r-integration in Eq. ([4]) we Fourier 
transform the Grassmann variables ipa using 

with the fermionic Matsubara frequency w„ = (2n + 
1)tt/(3 and obtain for the effective average action 



1 



2 ^ 

n—~oo 



*J,(ic^„E + C/f^)«'„. (9) 



Here E denotes the unit matrix. Moreover, we assume a 
regulator which is diagonal in Matsubara frequency space 



E 



i^*ti?,^„E'vI/„. 



(10) 



Inserting Ffc and Rk into Eq. ([T]) and comparing terms 
up to order zero in ip we obtain the following equation 

d_ 



2/3. 



(11) 



with K = {iujn + Uj^ )E' + i?fc.„E. From this equation 
we extract a set of ordinary differential equations for the 
various coupling constants by comparing the coefficients 
of the various monomials of Grassmann variables. In 
order to do this, it is necessary to invert the Grassmann 
matrix K. The matrix is inverted using the formula 



= {Ko + KiY 



N 

E(- 

m— 



-l\7n 



(12) 



where Kq is chosen to be the body of K (Grassmann 
scalar part of K) and Ki its soul (Grassmann non-scalar 
part). That the sum in Eq. only runs up to N is due 
to the fact that each matrix element of Ki is bilinear 
in the Grassmann variables, so that after TV factors of 
Ki the sum must terminate. Obviously, the inverse of 
K does not exist if the determinant of the body of K 
vanishes. In this case, the renormalization scheme is not 
well defined. 

At T = the sum over the Matsubara frequencies in 
Eq. ITTI) converts into an integral 



-y 



2^' 



(13) 



This integral is most easily evaluated for the sharp cut-off 
regulator 



(14) 



with C a suitably chosen large constant. This regulator 
fulfills the general requirements that it vanishes for k = 
and it dominates the effective average action for fc — s- 00. 
Furthermore, it has the advantage that the integration 
over the Matsubara frequencies can be done analytically. 
Some technical issues related to this integration will be 
discussed in Appendix A. 

An often used approach to the fermionic functional 
renormalization group applies a hard cutoff in such a 
way that the free propagator at large k is suppressed, 
and propagation is gradually 'switched on' when k is re- 
duced (see e.g. Refs.lTolandlsh. In the following we shall 
see that this (standard) method and the cutoff employed 
here yield identical sets of flow equations. 

With the results from Appendix A, we can write 
Eq. (fTTj) for the hard cut-off regulator (fTl]) in the fol- 
lowing form 



dk""''' - At: 



h E logdetfiG,T^(^A) 



X=±k 

N 



(15) 



dk 



(f^f-oo,..) - -^E E ^Tl■{~-MkGki^^)Y 



Att ^ — ' ^ — ' m 

7n—l X—±k 



(16) 



Eq. (|16p represents a finite set of coupled ordinary dif- 
ferential equations for the running couplings. The block- 
diagonal 2A^ X 2A^ matrix Gfe(iA) is given by 



Gk{iX) 











bkitX) 
bk{-iX) 



(17) 

and the matrix elements of the N x N submatrix bk{iX) 
are essentially determined by the self-energies 



(18) 



The 2N X 2N matrix Mk represents the soul of J/^^'e' and 
contains all running couplings except the self-energies. 

Eqs. (fT5|) and ([TB)l are the main results of this section. 
From these results, the flow equations for the various 
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running couplings can be directly extracted up to the de- 
sired order by comparing the coefficients of the different 
monomials of Grassmann variables. As will be discussed 
further below, up to second order in to, Eqs. (jlSp and 
(|16p are equivalent to results obtained in the literature 
using a hard cutoff on the free propagator. The advan- 
tage of Eqs. and (fTB|) is that they allow an easy and 
systematic construction of higher order contributions. 

In order to illustrate the use of Eqs. and 
we extract the flow equations for N = 2 assuming the 
symmetry ai2 — 021 and purely real running couplings. 
In this case, the submatrix gk{i^) is given by 



1 / 622 -bi2 

611&22 - 6?9 V ~^i2 ^11 



(19) 



with bap defined in Eq. ()18|) . The matrix Mk can be 
directly read off from Eq. ([7]) 



/ 



M, 



k 



ai2i2,fc 



*O10O2 





v 









001002 
002 002 
^02 001 



f01<P02 





/J02W1 
-001001 

- - - - /,Sl0O2 

-0Si0oi 
(20) 

After some algebra one obtains the following flow equa- 
tions for the various running couplings 

1 



= -Rclog (611622 -fo?2) A' 



ai2i2 



TT 



'1212 



Re- 



hiib 



'1212 



Re- 



22 - 
1 



6?2' 



(21) 



^1212 



611622 

6^1622 + 



&?2 



6116^2 - 26126*2 



27r (61162 



&?2) (&n&22 



The prime indicates a derivative with respect to fc, and 
we do not indicate explicitly the k dependence of the 
running couplings. 

For N = 2 the result Eq. pT|) is complete within the 
static approximation. This means, that due to the fi- 
nite number of Grassmann monomials which can be con- 
structed from four Grassmann variables, there will be 
no higher order equations. The right hand side of the 
equations for arises from m = 1 in the sum over to 
in Eq. p^ . while the right hand side of the equation 
for ai2i2 is the contribution from to = 2. Eqs. (|2ip are 
equivalent to Eqs. (26) and (27) in Ref.|3, where they are 
used to study electron transport through a single quan- 
tum dot. Moreover, the equations for in Eq. (|2ip are 
equivalent to Eq. (1) in Ref. 3 which was used to show 
that there are unusual correlation induced resonances in 
a polarized double dot system. (In order to show the 
equivalence one must set 01212, fc to its initial condition 
01212, fc = U.) 

We now discuss the structure of the set of differen- 
tial equations obtained from Eq. p6|) for larger N using 
A'' = 4 as an example. As can be easily established, the 



number of differential equations grows exponentially. In 
particular, for A" = 4 we find 35 coupled equations corre- 
sponding to the number of Grassmann monomials we can 
build from 8 Grassmann variables in a system with par- 
ticle conservation. If symmetries between these Grass- 
mann monomials are considered the number of equations 
reduces, e.g. for a system of four electrons with spin sym- 
metric interactions we only obtain 26 independent equa- 
tions. It turns out that the resulting (complete) set of 
equations is quite formidable and can only be generated 
with the help of computer assisted symbolic computation. 
Nevertheless, some features of this set of equations can 
be established in general. These follow from the struc- 
ture of Eq. p6p : A term of order to in the sum over to 
in Eq. p6|) only contributes to equations for coefficients 
with at least 2to indices. Such a term contains to prop- 
agators Gk- E.g., for A^ = 4, the term ni ^ 4 in the sum 
over TO only contributes to the equation for the running 
coupling ai234i234. Of course, the equation for this par- 
ticular running coupling also receives contributions from 
TO = 1,2,3. 

To further illustrate the use of Eq. P^ . we derive equa- 
tions suitable for the description of a one dimensional 
chain of nearest neighbor coupled spinless electrons (see 
e.g., Ref. [13) ■ To this end we first need to write down 
the effective potential in terms of the running couplings 

N 

N 

+^X!'^i'^j'^i+i'^i+i (22) 

with periodic boundary conditions 0o — ipN- Here, for 
simplicity, it is assumed that the density-density interac- 
tion strength U does not renormalize and remains at its 
initial value. In order to obtain the flow equations for 
the self-energies ajji , we consider the term to = 1 in the 
sum over to in Eq. p6p . One easily finds 

TV 

TiMkGk [5j_ij„i(iA) +.gj+ij+i(iA)] 05j-0oj 

- [Sij+i(«'^)0oj0oj+i +5j,j-i(*^)0Sj0oi-i] 
+ {i\ ^ -iX). (23) 

Inserting this result into Eq. p6p and comparing corre- 
sponding terms yields the fiow equations 



\=±k 



^ Yj fjj±i(«A) 



2lT 



(24) 



\=±k 



This result agrees with the one obtained using a hard 
cutoff on the free propagator; it is used to study e.g. 
persistent currents in mesoscopic rings in Ref. 
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FIG. 1: (color online) Polarized double quantum dot coupled 
to left and right leads. The arrows between the dots and the 
leads symbolize the tunnel couplings t''j. The arrow between 
the dots symbolizes interdot Coulomb interactions. 



III. TRANSPORT THROUGH CORRELATED 
QUANTUM DOTS 

In the following we study transport through various 
quantum dot (QD) systems using the "effective poten- 
tial" (i.e. static) approximation introduced in the pre- 
vious section. We will emphasize the role of higher or- 
der contributions beyond m = 2 in the sum over m in 
Eq. (|16p . Up to TO — 2, the equations resulting from 
Eq. (fTB|) can be shown to be equivalent to the flow equa- 
tions employed in Ref. 0. In order to elucidate the role of 
higher order contributions, we will study similar param- 
eter sets as considered in Ref. In some cases, calcu- 
lations using Wilson's numerical renormalization group 
are available for comparison. 

The standard model Hamiltonian for the description of 
QD systems (depicted schematically in Fig. [T|) consists 
of three essential terms: the dot Hamiltonian _ffD, the 
Hamiltonian of the non-interacting electrons of the leads 
He and the interaction iJoE between the QD and the 
leads, 

H = Hu + He + Hoe- (25) 
The dot Hamiltonian Hd is given by three terms 

jo- 

+ 1 E UJf'njany.' (26) 

where dj^ (dja) creates (annihilates) an electron in state 
j with spin cr, respectively. The site occupancy operator 
is denoted by nj^r = dj^dj^r. The dot levels eja may be 
shifted by application of a gate voltage Vg. The matrix 
elements tjji represent hopping amplitudes between dot 
levels and U^^f describes electron-electron interactions. 
Of course, J7// = 0. 

The non-interacting lead electrons created (annihi- 
lated) by c'j^i (cjai) at the Fermi surface of the right 
lead {I = R) and left lead (Z = L) are described by the 



Hamiltonian 

He = -Th Y '^]ai(^U+i)<yi + ^-c- (27) 

with a hopping amplitude Th between neighboring sites. 
The coupling of the leads to the quantum dot levels is 
modeled by a tunneling Hamiltonian 

H^^^-Yt\cl^^d,,„ + h.c. (28) 

jcrl 

where electrons can tunnel from or to the dot levels from 
site of the right or left lead. The tunneling matrix 
elements t^- are assumed to be real in the present paper. 
They determine the broadening of the dot levels as will 
be discussed in more detail below. 

The Hamiltonian H is translated into an action suit- 
able for a coherent state path integral formulation (see 
e.g., Ref.fisl) essentially by replacing creation and annihi- 
lation operators by corresponding Grassmann variables. 
Physical observables at temperature T = are obtained 
from the renormalized Green's functions of the correlated 
dot system (see Eq. ([T7|) ). 

gi^{iuj) = + aij,k=o- (29) 

In a static approximation, the aij^k=o do not depend on uj. 
The dot Green's functions g are N x N matrices where 
is the total number of electronic states included in 
the description of the quantum dot system. The matrix 
elements of the dot Hamiltonian eja + Vg, tjj> , and C/^y' 
provide the necessary initial values for the corresponding 
flow equations. All other flow equations start the flow 
from an initial value of zero. 

Since the dot system is coupled to the leads via the 
tunneling Hamiltonian -ffoE, the dot levels acquire an 
imaginary part which is determined using a standard pro- 
jection of the full Hamiltonian on the dot Hamiltonian as 
described in detail e.g. in Refs.[3 andllSl. 

gij^{iuj) = iujSij + ay,fc=o - isgn(u;)A.y (30) 

with the imaginary part given by 

A., = TTpo E (31) 

l=L,R 

in terms of the density of states of the lead electrons at 
the Fermi level po ~ l/(7'"''h) (in the wide band approx- 
imation). Usually, the absolute values of the tunneling 
matrix elements are given in terms of the parameters 

From this propagator, the average dot occupancies 
(rija-) at r = may be easily calculated using 

M = TT / due'^^^giiLu). (32) 

where the exponential factor enforces convergence of the 
integral at large uj. To relate the dot propagator to the 
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conductance is quite an intricate problem which has been 
addressed using both the Landauer-Biittiker— and the 
Kubo formahsni^^.. At T = 0, one obtains for the con- 
ductance Q the rather intuitive result 



(33) 



where the sum runs over all dot levels connected to the 
leads via non- vanishing tunneling matrix elements i^. 
The phase of the sum in Eq. (j33|) is known as the trans- 
mission phase. This quantity is accessible experimentally 
and will be also briefly discussed in the following. 



the figure clearly shows its importance (see also Ref. UM ) ■ 
We would like to emphasize, that within the static ap- 
proximation, there are no higher order vertex functions 
to be considered, since the set of equations (pij) is com- 
plete. It is therefore possible to apply these equations at 
even larger Coulomb interactions. One finds that with 
increasing Coulomb interactions the two outer Coulomb 
blockade peaks seen in Fig. [2] shift further outwards and 
the exponentially sharp correlation induced peaks in the 
center shift exponentially close to zero. The dot level 
occupancies displayed in Fig. [2] show the characteristic 
charge oscillations discussed in detail in Refs. Handi. 



B. Side coupled double dot 



A. Polarized double dot 

We first review calculations for a polarized double dot 
corresponding to two states with eio- = £20- = 0. Each 
state is occupied by only one electron as depicted in 
Fig. [1] The system is spin polarized, which may be 
achieved by putting it into a strong magnetic field. Stud- 
ies of this system using different methods were presented 
e.g., by Konig and Gefen^, Sindel et al^ as well as Meden 
and Marquard^. 

The effective potential for this problem corresponds to 
Eq. ([3]) and the flow equations for the running couplings 
are given by Eq. (PTjl . It is very easy to solve these equa- 
tions numerically. In Fig. [2] we show results for the pa- 
rameter set given in the figure caption. This parameter 
set with the rather large Coulomb interaction U/T = 8 
and ti2 — has been taken from Ref.fsl. (For convenience, 
here and in the following all energies will be given in units 
of r = r'.) The figure shows Coulomb correlation 
induced resonances close to Vg =0. The dashed curves 
agree with the results of Ref. H, which were calculated 
neglecting the renormalization of the Coulomb interac- 
tions. The full curves include this renormalization, and 



Now we consider a side coupled double dot with spin as 
depicted in Fig.|31 Only one dot is coupled to the leads di- 
rectly. For this system, the set of coupled flow equations 
is much more complicated than for the polarized double 
dot briefly reviewed above, and higher order equations 
m > 2 play a role. In fact, assuming spin symmetric in- 
teractions, we have to solve a set of 26 coupled equations 
as was pointed out before. We will investigate two phys- 
ically different cases: large interdot hopping ^12 > F and 
small interdot hopping ti2 < F. In both cases we study 
the influence of intra-dot Coulomb interactions U , which 
are chosen to be equal on both quantum dots. Inter-dot 
Coulomb interactions are set initially to zero, but they 
evolve to finite values during the renormalization run. In 
fact, all 26 running couplings (vertex functions) consid- 
ered here renormalize to non-zero values in certain gate 
voltage ranges, some of them attaining extremely large 
values. Side coupled double dots have been investigated 
by Cornalia and Grempel^- and Zitko and Bonca— using 
the NRG as well as Karrasch, Enss and Meden^ using 
the fRG truncated at m = 2. 

We start with a discussion of dot systems with large in- 
terdot hopping described by the parameter set ii2/r = 4, 
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-10 10 -10 10 -10 10 

Vg/r Vg/r Vg/r 

FIG. 3: (color online) Side coupled double dot with large interdot hopping: ti2/r — 4, U/F — 8, Ti = Ff — F/2, full lines: 
truncation order m = 4, dashed lines: truncation order m — 2, (a) conductance Q, (h) occupancy of directly coupled dot ni, (c) 
occupancy of side coupled dot n2 . 




FIG. 4: (color online) Side-coupled quantum dot device. 

U /T = 8, and eio- — e2a — and other parameters given 
in the caption of Fig. [3l The same parameter set was 
discussed in Refs. i and ll9j. Rough estimates discussed 
in Ref. [T^ suggest that one should observe conductance 
peaks located at Vg « ±(C//2 + ^12) with a width given 
by U/2. This is indeed seen in Fig. [S^a). In this fig- 
ure, we compare calculations truncated at different flow 
orders. The dashed curves show calculations truncated 
at TO = 2 as in Ref. U. The full curves include all or- 
ders up to TO = 4, which is the maximum order possi- 
ble in this system. We observe that the inclusion of the 
higher orders has the effect of somewhat increasing the 
width of the peaks in accordance with the NRG calcula- 
tions presented in Ref. llH. The dot occupation numbers 
[Fig. [3I^b)and (c)] show slightly more pronounced shoul- 
ders if higher order effects m > 2 are included within the 
gate voltage ranges of peak conductance. 

Now we increase the intra-dot Coulomb interaction to 
C//F = 12. According to the discussion presented in 
Ref. 1^1, the peaks should become a more box-like in 
shape. This is indeed seen in Fig. OJa), but only if the 
higher order contributions m > 2 are included. A cal- 
culation truncating the set of equations at m = 2 shows 
conductance peaks with a Lorentzian shape in disagree- 
ment with NRG calculations^. The conductance peaks 
shown in Fig. [5l^a) are not completely flat as suggested 
by the NRG calculations. This probably points to effects 
beyond the static approximation employed here. The oc- 
cupancies of the two dots are shown in Fig.[5^b) and (c). 



Here the effect of higher order contributions is quite im- 
portant as well. The single occupancies rii even show 
a non-monotonic behavior in the gate voltage ranges of 
peak conductance. Similarly, the transmission phase de- 
velops a shoulder in these gate voltage ranges which is 
not seen in calculations truncated at to = 2. 

If one increases the intra-dot Coulomb interaction fur- 
ther, e.g., to U/T = 16, one observes conductance peaks 
which are flat, but within this flat region the renormaliza- 
tion procedure employed here becomes ill defined. Tech- 
nically, this can be traced to the fact that somewhere 
along the renormalization flow the renormalized self en- 
ergies of the indirectly coupled dot vanish. At this point 
the body of the matrix K which was introduced in the 
previous section vanishes, all flow equations develop a 
singularity, and the renormalization flow stops. The sin- 
gular renormalization flow likely indicates a failure of the 
static approximation. 

Let us now discuss a physical situation with small 
inter-dot hopping: Fig. [6] shows a calculation using the 
parameters ti2/T = .2 and U/T = 2, e\a = £20- = 0. This 
parameter set was also investigated in Refs. and [H. As 
was pointed out by the authors of Ref. 0, the calcula- 
tion truncated at to = 2 yields a conductance minimum, 
which is too flat in comparison to NRG calculationai^. 
As is obvious from Fig. inclusion of higher order ef- 
fects improves this situation. However, comparison with 
the NRG data presented in Ref. Il^ still shows quanti- 
tative differences in this gate voltage region. Since our 
calculations are not truncated within the constraints of 
the static approximation , these differences must be at- 
tributed to effects beyond the static approximation. 

Of course, the structure obtained for the conductance 
and the occupancies as a function of the gate voltage can 
be related to the calculated spectrum of the dot system. 
In Fig. [7] we show the energy level for a spin up electron 
on the directly coupled dot given by an and the energy 
level of a spin up electron on the side-coupled dot 033 
for the parameter set used in Fig. [S] Again we compare 
calculations truncated at to = 2 with a full calculation 
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FIG. 5: (color online) Side coupled double dot with large interdot hopping: ti2/T = 4, U/T = 12. Other parameters and notation 
as in Fig. [21 
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FIG. 6: (color online) Side coupled double dot with small interdot hopping fi2/r = .2, U/T — 2. Other parameters and notation 
as in Fig. [3 



in the static approximation. Obviously, truncating at 
m = 2 yields significantly too small energy eigenvalues in 
particular for the indirectly coupled dot. 

If we increase the intra-dot Coulomb interaction from 
U/r = 2 to t7/r = 2.5, we observe a picture similar to 
the one shown in Fig. [6] with the amplitude of the charge 
oscillations on the dots significantly increased. Already 
for U/r = 3, the set of flow equations develops singular 
behavior as was discussed above for the case with large 
inter-dot hopping. This again indicates that a truncation 
scheme based on a static approximation fails for large 
Coulomb interactions . 



C. Parallel coupled double dot 



ti2 = 0. This parameter set has also been investigated 
in Ref. 4. Due to the large Coulomb interaction, correla- 
tion induced resonances are observed similar to the spin- 
polarized system briefly reviewed in section Fill Al How- 
ever, as is seen in Fig. [HI calculations truncated at m = 2 
show significant quantitative differences in comparison 
to the full calculation up to m = 4: The resonances are 
much more pronounced and the valley between these res- 
onances is somewhat narrower. This trend continues, if 
we increase the Coulomb interactions to U^p /F = 4. Re- 
sults are shown in Fig. [TOl The amplitude of the charge 
oscillations of the dot with the stronger tunnel coupling 
to the leads is significantly underestimated in a calcula- 
tion neglecting higher order contributions. 



Finally we consider a parallel coupled double dot as 
depicted in Fig. [S] The set of equations to be solved is 
the same as for the side-coupled dot. 

We study a parameter set where all t\ are chosen pos- 
itive with the tunneling strengths F' given numerically 
in the caption of Fig. M We choose U^f'/T = 2 for afl 
intra-dot and inter-dot interaction matrix elements and 



Of course, the question arises, to what extend the 
static approximation is still valid in this parameter re- 
gion. In order to answer this question calculations be- 
yond the static approximation must be performed. Such 
calculations are presently under way. 
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FIG. 8; (color online) Parallel coupled quantum dot. 

IV. CONCLUSIONS 

The functional renormalization group has been pro- 
posed as a useful new tool for the investigation of in- 
teracting mesoscopic systems'^-^. The advantage of the 
method in comparison to NRG or DMRG calculations 
appears to be its rather moderate numerical cost. How- 
ever, this advantage is due to some rather drastic ap- 
proximations, which are not easily controlled. In this 
paper we carry out the first steps that are needed to go 
beyond the standard set of approximations usually em- 
ployed to solve the fermionic fRG in practice. Still within 
the constraints of a static approximation we develop a 
method which systematically generates the complete set 
of ordinary differential equations corresponding to the 
fRG. However, the method proposed here is not limited 
to the static approximation. It can be extended straight- 
forwardly to include non-static effects e.g. by including a 
wave function renormalization in the kinetic energy term 
of the action. 

We demonstrate by means of a number of calcula- 
tions of transport properties of quantum dot systems that 



inclusion of higher order equations definitely improves 
fRG calculations for large Coulomb interactions. This 
is judged by a comparison to NRG results. However, if 
one increases Coulomb interactions beyond some critical 
value, the static approximation fails as is signalled by sin- 
gularities developing during the renormalization flow. In 
order to overcome such problems and to further improve 
the agreement between NRG and fRG calculations it is 
necessary to go beyond the static approximation. Work 
in this direction is presently under way. 
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APPENDIX A: INTEGRAL OVER THE 
MATSUBARA FREQUENCIES 

In order to perform the rather singular integral over 
the Matsubara frequencies ui for the hard cutoff regulator 
Eq. ([M]) we need the formula^i 

6{x - y)f{9{x - y)) S{x - y) f dsf{s). (Al) 

Jo 

The w-integration leading from Eq. pTjl to Eq. ITSl) is 
then performed as follows 



Tr / — {G-\iuj) + cke{e ^ uj'^m ^ ^cke{k^ - oj^) 

^ ' ok 



2tt 



^+Tr r ^(G-\iiu) + Cke{k^~Lu^)Ey'2Ck' 

7_oo 27r 



S{k 



(A2) 



10 




v,/r 



-6 -4 -2 2 4 6 

v,/r 



-6 -4 -2 2 4 6 

v,/r 



FIG. 9: ("cotor onlme) Parallel coupled double dot ii2/r = 0, U/T = 2, Tf/F = .5, Tf /F = .25, F^- = 0.07, = 0.18. Oi/ier 
notation as in Fig. [21 
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FIG. 10; (color online) Parallel coupled double dot ti2/r = 0, U/F = 4, Ff/F = .5, Ff /F = .25, = 0.07, Ff = 0.18. Other 
notation as in Fig. [3 



■uS,,^-')E'. The quantity w^^"* rep- 



with G """(iw) = {iuj r iij, 

resents the body of C/^: . In the first term we aheady 
performed the hmit C oo. The second term is now 
evaluated with the help of Eq. (|A1|) 



Tr 



duj 

Ck {S{k-uj))+S{k + uj)) 

rC 



= -j^Yl Trlog^ic-'liA)] +0(logC) 



A=±fc 



Therefore we obtain the equation 



\=±k 



The infinite constant gets absorbed into a suitably mod- 
ified initial condition. 

With similar steps as above one obtains the result given 
in Eq. (fTB)) . Here it is important to use the fact that the 
factors under the trace may be permuted cyclically. 
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